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ABSTRACT 

Anisotropic gravitational radiation from a coalescing black hole binary is known 
to impart recoil velocities of up to ~ 1000 km s~^ to the remnant black hole. In this 
context, we study the motion of a recoiling black hole inside a galaxy modelled as 
an Hernquist sphere, and the signature that the hole imprints on the hot gas, using 
N-body/SPH simulations. Ejection of the black hole results in a sudden expansion of 
the gas ending with the formation of a gaseous core, similarly to what is seen for the 
stars. A cometary tail of particles bound to the black hole is initially released along 
its trail. As the black hole moves on a return orbit, a nearly spherical swarm of hot 
gaseous particles forms at every apocentre: this feature can live up to ^ 10^ yr. If 
the recoil velocity exceeds the sound speed initially, the black hole shocks the gas in 
the form of a Mach cone in density near each super-sonic pericentric passage. We find 
that the X-ray fingerprint of a recoiling black hole can be detected in Chandra X-ray 
maps out to a distance of Virgo. For exceptionally massive black holes the Mach cone 
and the wakes could be observed out to a few hundred of Mpc. Detection of the Mach 
cone is found to become of twofold importance: i) as a probe of high-velocity recoils 
and ii) as an assessment of the scatter of the Mbh — ^buige relation at large black 
hole masses. 
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1 INTRODUCTION 

Today, black holes (BHs) with masses in excess of 10^ 
appear to be ubiquitous in bright galaxies (Richstone et al. 
1998; Decarli et al. 2007) and scaling relations between the 
BH mass and the underlying galaxy indicate unambiguously 
that BHs evolve in symbiosis with their hosts (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000; Hopkins et al. 2006). 

According to the current paradigm of structure forma- 
tion, galaxies often interact and merge as their dark matter 
halos assemble in a hierarchical fashion, and the BHs, incor- 
porated through mergers, are expected to grow, evolve and 
pair with other BHs ( Volonteri et al. 2003) . The formation 
of super- massive BH pairs thus appear to be an inevitable 
and natural consequence of galaxy assembly. 

In our local universe, the dual radio source 3C 75, at the 
centre of Abell 400, is a clear, albeit rare example of BH pair- 
ing as it hosts two massive BHs displaying prominent radio 
jets (Owen et al. 1985). Chandra observations of 3C 75 have 
revealed the occurrence of two active nuclei at a projected 
separation of 7 kpc (Hudson et al. 2006) providing first ev- 



idence of their current coupling. An even more remarkable, 
and still unique example is the case of the elliptical galaxy 
0402+369 (Rodriguez et al. 2006). VLBI observations high- 
lighted the presence of two compact variable, fiat-spectrum 
radio sources at a projected separation of only 7 pc, suggest- 
ing that the two massive BHs that power the radio emission 
form a binary. These two cases indicate that BH pairing oc- 
curs and proceeds from the large-scale of a merger (up to 100 
kpc) down to the scale where the BHs form a binary in close 
Keplerian relative orbit (of a few pc). In the interaction with 
the stars (Milosavljevic & Merritt 2001; Yu 2002; Berczik et 
al. 2006; Merritt & Milosavljevic 2005; Sesana et al. 2007) 
and/or gas (Armitage & Natarajan 2002; Escala et al. 2004; 
Escala et al. 2005; Dotti et al. 2006; Dotti et al. 2007; Mayer 
et al. 2007) that surround the binary, the BHs lose orbital 
angular momentum, and if the process continues down to a 
scale of a few milliparsec, gravitational wave emission drives 
the BH in-spiral down to coalescence, causing the formation 
of an heavier remnant BH. 

Recently, a major breakthrough in numerical relativity 
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has allowed to trace for the first time the BH binary evolu- 
tion down to coalescence (Pretorius 2007), in the strong field 
regime imposed by general relativity under arbitrary condi- 
tions (Campanelh et al. 2007A; Campanelh et al. 2007B; 
Bruegmann et al. 2007; Baker et al. 2007; Baker et al. 2008; 
Herrmann et al. 2007; Herrmann et al. 2007; Herrmann et 
al. 2007; Koppitz et al. 2007). These studies have revealed 
that spinning BHs emit an anisotropic beam of gravita- 
tional radiation, and in response to this asymmetry the rem- 
nant BH receives a recoil (due to linear momentum conser- 
vation) that can displace it from the central parts of the 
galaxy. In more detail, BHs with similar masses and spin 
vectors in non-generic alignments relative to the orbit pro- 
duce anisotropic patterns of radiation via spin-orbit coupling 
that lead to recoil velocities Vrec in the range of < 200 km s~^ 
to 2000 km s~^, and for particular spin-orbit configurations 
to a recoil as high as 4000 km s~^ (Campanelh et al. 2007A; 
Campanelh et al. 2007B; Baker et al. 2008; Schnittman & 
Buonanno 2007). These "natal kicks" are in the range of 
the typical dispersion/rotation velocities of galaxies or even 
higher, so that the remnant BH can either be significantly 
displaced away from the galactic centre where the merger 
occurred, or be ejected when Vrec exceeds the escape speed 
from the galaxy (typically of ~ 1000 km s~^ for massive 
galaxies) . 

Both a BH that remains or a BH that escapes can leave 
a "sign" on the underlying galaxy with observable conse- 
quences. As an example, a BH ejected from its host galaxy 
can carry a punctuated accretion disc that lights the BH on, 
as an active X-ray source, eventually deprived of its under- 
lying galaxy (Loeb 2007). Can a BH retained in the galaxy 
carry other observable features? This is the question that 
we address in this paper. 

A kicked BH, retained inside the host galaxy, moves 
along a radial elongated orbit; it both explores the galaxy 
periphery and cross the centre a few times, periodically. Dy- 
namical friction against stars and gas causes its radial orbit 
to decay over a time-scale that depends on the recoil velocity 
and underlying background density. 

The BH damped oscillatory motion, in a pure stellar 
background, has been investigated by a number of authors 
(Merritt et al. 2004; Boylan-Kolchin et al. 2004; Gualandris 
& Merritt 2007) . After an early phase of braking under the 
action of dynamical friction and a number of recursive pas- 
sages across the centre of the galaxy, the orbital decay cause 
the motion of the BH to be confined into its gravitational 
influence radius. From this point on the BH starts oscillat- 
ing together with the core around the common centre of 
mass and it eventually reaches thermal equilibrium with the 
stars when the oscillations decays down to Brownian level 
(Gualandris & Merritt 2007). 

The interaction of the massive BH with stars has no neg- 
ligible effects on the underlying equilibrium (Boylan-Kolchin 
et al. 2004; Merritt et al. 2004; Gualandris & Merritt 2007), 
the most interesting being the formation of a stellar core. As 
pointed out by Boylan-Kolchin et al. (2004), the instanta- 
neous removal of the BH from the galaxy centre, at the time 
of its ejection, causes a sudden decrease in the gravitational 
potential so that stars driven away from virial equilibrium 
respond dynamically expanding. In addition, when the BH 
returns to the centre, it transfer orbital kinetic energy into 



the background under the action of dynamical friction, and 
this energy deposition causes further stellar expansion. 

In this paper we address a complementary problem, i.e. 
the effect that a recoiling BH has on the hot gaseous com- 
ponent of an early type galaxy host. Local temperature and 
density perturbations excited by the super-massive BH while 
it is travelling across the galaxy are expected to lead to ob- 
servable changes of the bremsstrahlung emission. We will 
specifically address the following issues i) how the struc- 
tural properties of the gas in the galaxy are modified by the 
ejection of the BH and its orbital decay; ii) which signa- 
tures the BH motion imprints on the gas; iii) the prospects 
for detection of these signatures by X-ray telescope such as 
Chandra. 

For this study, we perform a suite of 8 simulations of gas 
rich elliptical galaxies, where BHs are ejected with different 
recoil velocities. The BH mass Mbh scales with the host 
mass according to the observed correlation (Haring & Rix 
2004). Since this correlation has an intrinsic scatter of at 
least 0.3 dex (see Lauer et al. 2007, Tundo et al. 2007 for 
a discussion of the scatter of the correlations at large BH 
masses), we also study the case of an over-massive BH with 
a mass a factor of 3 larger than predicted by the best fit 
value of the Mbh — Mbuige relation in Haring & Rix (2004). 

The outUne of the paper is as follows. In Section 2 we 
introduce the physical model used in our simulations. In 
Section 3 we describe the changes in time of the density and 
temperature induced by the ejected BH, inferred consider- 
ing different recoil velocities, and for a BH mass following 
the Mbh — Mbuige relation and its scatter. In Section 4, we 
explore the BH detectability with Chandra. In Section 5 we 
discuss our results. 



2 SIMULATIONS 

We perform our simulations using the N-Body/SPH code 
GADGET (Springel et al. 2001). The elliptical galaxies are 
modelled as an Hernquist sphere (Hernquist 1990) with scale 
radius a and total mass Mtot- The mass density, cumulative 
mass profile and escape velocity are given by: 



p{r) 



Mtc 



27r r(r -h a)^ 



M(r) = Mtot 



(r -h a)2 



Vesc = (2GMtot/a)'/'. 



(1) 

(2) 
(3) 



The total mass of the galaxy is Mtot = 10^^ M©, with 
a gas fraction /gas = Mgas/Mtot =0.1; the gas follows the 
same density profile of the stellar component. We explore 
two different scale radii: a equal to 4 and 8 kpc , and two 
different BH masses, 2 and 6x10^ M© (see Table 1). The 
lightest BH follows the Mbh — Mbuige relation as in Haring 
& Rix (2004), while the heaviest explores the scatter at high 
BH masses (Lauer et al. 2007, Tundo et al. 2007) where we 
expect the imprints of the recoils to be the largest. 

For the colUsionless component we use 10^ particles 
with a gravitational softening of 100 pc to prevent numeri- 
cal collisional relaxation. In the runs with the light (heavy) 
BH we use 2 x 10^ (10^) gaseous particles with softening 
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hgas = 10 (50) pc; the number of gas particles used to aver- 
age the hydrodynamical properties is 40, so that our mass 
resolution is ^ 2 x lO^M© (10^ M©). Each gaseous particle 
is evolved along a poly trope P oc with 7 = 5/3, corre- 
sponding to adiabatic evolution of an ideal gas. Compres- 
sional work and shock heating are included in the energy 
equation. The gravitational softening of the BH equals that 
of the gas, /ibh = /igas- We ensure that the dynamical in- 
fluence of the hole on the gas is resolved, by checking that 
hBU <^ Hnf — a(2 MBH/Mtot)"*^^^- Here the influence radius, 
Tinf is defined as the distance from the galaxy centre that 
encloses a total mass equal to 2 Mbh • For the less massive 
BH: 

" [j^c ) V 2xlO^Mo M.ot ) 

We set the initial condition for the BH plus galaxy sys- 
tem following the prescription by Hernquist (1993), and we 
evolve the equilibrium model for 1 Gyr to verify its stability. 
After the ejection, we halt our simulations when the BH has 
settled into its Brownian motion. 

We explore the ejection of the BH for different recoil 
velocities Vrec ranging from 0.2 to 0.95 the escape velocity, 
Vesc- In the simulations with a = 4 (8) kpc the escape ve- 
locity at r = is 1490 (1054) km s~^ and the central initial 
sound speed is 400 (300) km s~^. The resulting initial Mach 
number, defined as = t'rec/cs,o ranges between 1 and 3.3. 
Table 1 lists physical parameters for all 8 models: a, through 
H. We will identify with lower-case (upper-case) letters the 
runs referring to the light (heavy) black hole. 
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Figure 1. Radial density profiles of gas, in units of the initial 
central density at the softening radius, as a function of the dis- 
tance from the centre of the galaxy in units of a for run b at 
different times. Solid line refers to the initial condition, dotted 
and dashed lines to the time at which the BH is at its first apo- 
centre and first pericentre respectively. Dot-dashed line refers to 
the final configuration when the BH orbit has decayed. Vertical 
line indicates the position of Hnf. Red dots refers to the run g 
where the BH have been instantaneously removed. 



3 BH INTERACTION WITH THE HOT GAS 

The recoiling BH reaches its first apocentre rapo ^ 
a [(1 — 'i^rec/'^esc)"^ " l] • Followlug the kick, the motion of 
the BH exhibits three distinct phases: (i) a decaying os- 
cillatory radial motion guided by dynamical friction; (ii) a 
phase characterised by BH oscillations well inside the galaxy 
scale radius followed by (iii) Brownian motion. These phases 
agree with what Gualandris & Merritt (2008) found for the 
collisionless case. 

In the following we investigate the direct influence that 
the ejection of the BH has on gas density and tempera- 
ture, looking at the evolution of both radial profiles and 
2-dimensional maps. 

3.1 Density and temperature radial profiles 

We now turn our attention to the global disturbance that 
the BH imprints on the gas. We first focus on a simu- 
lation with Mbh = 2 x 10^ M© (run h, M = 3), and 
analyse how the density and temperature profiles change 
during the entire run. In run b the orbit decays into the 
Brownian regime after ~ 458 Myr (i.e., ~ 20tdyn where 
tdyn = (3^2aV(GMtot)'/'). 

Figure 1 shows the evolution of the radial density pro- 
file. When the BH is at its first apo-centre (22 Myr after the 
ejection) the central density drops creating a core within 
^ 2rinf. The density keeps decreasing until the BH sets 



into its Brownian motion. The final central gas density is 
~ 9 — 10 times lower than the initial value, and the final 
gaseous core has a size ^ 4rinf. 

The evolution of the radial density profile is caused by 
the same two processes outlined by Boylan-Kolchin et al. 
(2004) and Merritt et al. (2004) for the stellar component. 
At first the density drop is due to the sudden removal of 
the BH gravitational field and not to dynamical friction: in 
response to this change, the stellar and gaseous components 
readjust to a dynamical equilibrium state with no BH. To 
check this hypothesis we carried out a test simulation where 
the BH was removed instantaneously from the galaxy (run 
g for Mbh = 2 x 10^ M©) and we let the system evolve for 
22 Myr. Density profiles for run b and run g at this time 
are indeed very similar, confirming this initial hypothesis. 
At later times, the density evolution is driven by dynamical 
friction, which injects energy into the background at each 
supersonic pericentric passage. When the orbital decay is 
completed, the combined action of the two processes ends 
with an erosion extending all the way to 4nnf. 

Figure 2 shows the evolution of the temperature profile 
in run b. The central temperature T decreases in response 
to the ejection of the BH as well. However, when the BH 
settles down in its Brownian motion, T returns to its virial 
value. 

The final gas density profiles normalised to the initial 
profiles are shown in Figure 3 for cases from a through F. 
Core formation occurs in all cases. For a fixed Mbh, the 
importance of the core increases with Vrec^ while for a fixed 
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Run 


Mbh 


'^esc 


'^rec 


M 


Mdef,star/MBH 


Mdef,gas/MBH 


Mdef,tot/MBH 


''"core/ CL 


Mgas,b/MBH 


a 


2 


1490 


500 


1.3 


1 


0.24 


1.24 


0.16 


0.13 


b 


2 


1490 


1200 


3 


3.2 


0.63 


3.83 


0.35 


0.01 


C 


6 


1054 


400 


1.3 


1.8 


0.19 


2. 


0.16 


0.3 


D 


6 


1054 


700 


2.3 


3.5 


0.35 


3.85 


0.275 


0.1 


E 


6 


1054 


900 


3 


3.8 


0.43 


4.2 


0.36 


0.061 


F 


6 


1054 


1000 


3.3 


4. 


0.53 


4.5 


0.41 


0.05 


g 


2 


1490 


oo 


oo 


0.73 


0.06 


0.79 


0.075 




H 


6 


1054 


oo 


oo 


1.53 


0.21 


1.75 


0.19 





Table 1. List of the simulations: label of the run; BH mass in unit of 10^ M©; escape velocity at r=0 and initial recoiling velocity in 
kms"-*^; initial Mach number A4; mass deficit in stars, gas and their sum in units of Mbh; core radius in unit of scale radius, a; mass of 
particles initially bound to the BH Mgas,b i^i units of Mbh- Run g and H correspond to the case of instantaneous removal of the BH. In 
run a, b, g (C, D, E, F, H) the scale radius is equal to 4 (8) kpc and the initial sound speed is 400 (300) kms"-*^. 




Figure 2. Radial temperature profiles of gas as a function of the 
distance from the centre of the galaxy in units of a for run b at dif- 
ferent times. Solid line refers to the initial condition, dotted and 
dashed lines refer to the first apocentre and pericentre, respec- 
tively. Dot-dashed line refers to the final configuration. Vertical 
line indicates the position of rinf. 



Figure 3. Final gas density profiles pgas as function of r/a, nor- 
malised to the density profile at t=0 Pgas,o for the two sets of run 
with different BH masses and different recoil velocities Vrec (in 
units of km s"-*^). Vertical lines refer to the influence radii of the 
two BH masses: light (heavy) curve for Mbh =2(6) xlO^M© 



Vrec it increases for more massive BHs, i.e., with the kinetic 
energy deposited by the BH. For low recoil velocities (M ^ 
1), the kinetic energy carried by the BH and deposited in its 
transit across the centre of the galaxy is low, and the major 
perturbation acting on the gas is related to the temporary 
removal of the BH at the time of coalescence. 

We compute final core radii, and mass deficits both for 
the gaseous (Mdef,gas) and stellar (Mdef,star) components 
(see Table 1). Mdef,star is always larger than Mdef,gas, be- 
cause in all our simulations the stellar mass is larger com- 
pared to the gaseous one. The total, star+gas, deficits span 
from 1.5 up to 4 Mbh, in agreement with the results from 
coUisionless simulations (Boylan-Kolchin et al. 2004; Gua- 
landris k Merritt 2007; Merritt et al. 2004). 



3.2 Density and temperature 2-D maps 

In this section we study the shape and extent of the density 
and temperature perturbations that the BH excites along 
its trail. We first discuss models with ^ 1, i.e. run a and 
C. During the orbital decay the BH is always surrounded 
by a hot spherical over-density of gaseous particles. The size 
of the spherical overdensity corresponds roughly to rh = 
GMbu/v^qc^ i.e. the radius within which particles remain 
bound after the BH is kicked. The density enhancement is 
illustrated in the upper panels of Figure 4 for both runs 
when the BH is at its first apocentre. 

Another characteristic feature in the density maps of 
these early dynamical phases, is a stream of gas particles lag- 
ging behind the BH trail. We argue that this gas is stripped 
material initially residing inside rinf. To test this hypothesis 
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Figure 4. Upper panel: color-coded gaseous density maps in loga- 
rithmic scale for run a (left) and C (right) (-Urec = 500 km s"-*^ and 
400 km s"-*^, respectively) when the BH is at its first apocentre. 
Brighter colors refer to higher densities that range between 0.1 
and 10~^ cm~^. Lower panel: same maps but where the gaseous 
particles initially bound to the BH are marked in green. Boxes 
are 2 a on a side (i.e. 8 and 16 kpc for run a and C, respectively). 
The BH initial orbit is on the x-axis and maps are projected in 
the x-y plane 



we mark, at the start of runs a and C, the particles gravita- 
tionally bound to the BH^, and then follow their dynamics. 
Figure 4 shows the density map of the gas surrounding the 
BH at its first apocentre. In the lower panels the particles 
originally bound to the BH are marked in green. Here a clear 
spatial coincidence between the stream and the marked par- 
ticles appears, thus supporting our hypothesis. 

As we move to higher velocities {A4 > 2), the mass 
of particles initially bound to the recoiling BH (Mgas,b) de- 
creases considerably as reported in Table 1. Again, this was 
previously noted for the stellar case by Gualandris & Merritt 
(2008). The quasi spherical swarm of hot gaseous particles 
around the BH becomes less relevant, and for the case with 
the highest speed (runs b and F) , it is no longer visible even 
soon after the ejection. The lack of particles around the BH 
for large Mach numbers is a consequence of our resolution 
limits: for kick velocities ^ 800kms~^ (700kms~^) and for 
the light (heavy) BH, the radius rb drops below our resolu- 
tion limit so that we are unable to resolve the initially bound 
particles anymore. 

New features appear when A4 > 2 initially: a steep, con- 
ical over density develops, i.e. a Mach cone becomes visible 
in the density maps, as shown in Figure 5. The opening angle 
of the cone is inversely proportional to Ai , and varies along 
the BH trajectory as both its speed and background tem- 
perature vary. The shape of the cone is sharpest at the first 
pericentre passage, where the contrast between the speed of 

^ We define particles gravitationally bound to the BH as done 
in Gualandris &; Merritt (2008): we calculate the relative energy 
between the BH and each particle. When this energy is negative 
we define the particle as bound to the BH. 



Figure 5. Upper panel: color-coded gaseous density maps in lin- 
ear scale for run b (left) and E (right) with J\4 = 3 when the BH 
is at its first pericentre. Density range between 0.001 and 0.05 
cm.~^q. Lower panel: color-coded temperature map in logarith- 
mic scale between 5 x 10^ K and 1.5x10^ K. Boxes are 4 and 20 
kpc on a side for run b and E, respectively. The BH initial orbit 
is on the x-axis and maps are projected in the x-y plane 



the perturber and the temperature of the gas is the highest 
(see Figure 2). The cone then weakens during the rest of the 
simulation as dynamical friction slows down the BH motion, 
and the temperature returns to its virial value. The Mach 
cone is clearly visible along the direction orthogonal to the 
BH velocity vector, and it can still be recognisable as long 
as the angle between the line of sight and the BH is below 
^ 45 degrees. 

All these features are transient and have a characteris- 
tic lifetime that depends on the BH-galaxy mass ratio, on 
the extent of the recoil velocity imparted at the moment 
of BH coalescence, and on the background. For low Mach 
numbers, M ^ 1, the over-density around the BH is more 
easily visible when the BH is at apocentre as the density and 
temperature contrast with the background are the highest. 
As the BH spends most of the time at apocentre, this over- 
density can be observed for 10 — 100 Myr. For high Mach 
numbers, i.e. recoil velocities near escape, the Mach cone 
becomes the main feature and it is seen during the first su- 
personic pericentres close to the central region of the host 
galaxy, within a scale ranging between 0.5-2 kpc. The cor- 
responding lifetime is a few tens of Myr. 

3.3 Gas fraction 

We now briefly discuss how the density and temperature 
perturbations excited by the BH travelling across the hot 
gas of the galaxy depend on the fraction of gas /gas relative 
to the stellar mass. Although we have modelled gas rich 
early type galaxies with /gas = 0.1, the dynamical evolution 
of the BH is still dominated by the stellar component. As 
a consequence, a further reduction of /gas does not alter 
the BH orbit or the characteristic decay time. As the gas 
is in virial equilibrium, its sound speed and temperature 



6 Devecchi, et al 



are determined by the overall gravitational potential, and 
thus do not vary for smaller /gas- Since the shape of the 
perturbation depends on the value of the Mach number and 
on the size of the BH influence radius, galaxy models with 
different gas fractions but same M and Mbh, differ only in 
the total X-ray luminosity (discussed in Section 4). We run 
a test simulation with /gas — 0.01, M — 2.3 and Mbh — 
6 X 10^ M0 to verify these scaling and find good agreement 
with our expectations. 



4 DETECTABILITY 

The motion of the BH perturbs the density and temperature 
distribution. Bremsstrahlung emission from the perturbed 
gas is expected to create X-ray features above the diffuse 
emission from the underlying hot gas in the galaxy. In this 
Section, we quantify the detect ability of signatures related 
to recoiUng BHs in X-rays. We create synthetic Chandra 
observations using the code X-ray MAp Simulator (Gardini 
et al. 2004; Rasia et al. 2008, X-MAS). The mock event 
files are generated in the soft band, [0.7, 2] keV, assuming 
the Response Matrix File and the Ancillary Response File 
of the Chandra detector ACIS-S 3. The images are created 
assuming an exposure time of 10 ks for the heavy BH and 
100 ks for the light BH. The galaxy-BH system is considered 
to be at the distance of Virgo (z=0.004), meaning that the 
field of view of our detector (8.3 arcmin) corresponds to 42 
kpc. The orbital plane of the BH is perpendicular to the 
X-ray image and oriented along the horizontal axis unless 
otherwise noted. 

The drop in central density is less dramatic for the light 
BH and the swarm of bound particles less evident. The X- 
ray contrast of the features produced over the underlying 
galaxy is consequently weaker than in runs with a heavy 
BH. Identification of recoil signatures on raw surface bright- 
ness maps is hard. To better highlight the influence of the 
BH in Chandra maps and to quantify its relevance over the 
background galaxy, we create "mirror images" . We rotate a 
galaxy image by 180 degrees around the axis orthogonal to 
the BH motion and then divide the original image by the ro- 
tated one (this technique is reminiscent of the "Asymmetry" 
index used in optical astronomy, e.g., Conselice 2003). 

Selected mirror images for the Ught BH are shown in 
Figure 6. When M — 1 (run a), the perturbation of the 
gas results in an X-ray enhancement off-set relative to the 
galaxy centre. Its concentration and spherical shape can be 
associated to the feature outlined in Figure 4. The three 
images for run b represent three different moments of the 
BH orbital decay. During the first oscillations the BH has an 
orbit with large radii (~ 8 kpc ). The residues of its passage 
are visible, as broad wings, up to 5 kpc from the centre. As 
the decay proceeds the orbital radius decreases and in this 
particular simulation the BH starts to precess after ~ 200 
Myr. The second mirror image is taken after 266 Myr when 
the motion axis is tilted by almost 45 degrees respect to its 
original orientation. Here the BH motion is traced by the 
two yellow spots oriented along the BH trajectory. The last 
image refers to the last phase of the BH motion when the BH 
moves in a short orbit across the centre and passes frequently 
through the same excited region. The X-ray enhancement 
increases the background emission of 30%. 




Figure 6. Mirror images for the light BH smoothed with a gaus- 
sian filter. Left-Top panel: Run a, when the BH is in its first apoc- 
entre. Right-Top, Left-Bottom and Right-Bottom panel: Mirror 
images for run b at three different moments of the BH motion 
(time after the ejection are shown in each panel). The color scale 
is linear between 0.8 and 1.2. 
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Figure 7. Simulated color coded maps of the X-ray emission de- 
tectable by Chandra (ACIS-S3 chip), for run D. The photon im- 
ages are computed in the [0.7 2] keV band, background subtracted 
and 1 arcsec binned. The angular size is 2.7 arcmin correspond- 
ing approximately at 14 kpc at the source redshift (z=0.004). The 
scale is logarithmic and the isophotal contours (shown in black) 
are logarithmically equi-spaced between 100 and 50000. Different 
panels correspond to first apocentric phases (denoted by Apo) and 
pericentric ones (denoted by Peri). The last pericentre showed is 
after 145 Myr, while the last apocentre is at 120 Myr. 
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Figure 8. Upper panels: Simulated color coded maps of the X- 
ray emission detectable by Chandra (ACIS-S3 chip), for run E 
for different angles (reported in each image) between the plane 
orthogonal to the line of sight and the BH motion. Lower panels: 
Mirror images for run E (with the heavy BH) smoothed with a 
gaussian filter for different rotation angles between the line of 
sight and the BH motion. The color scale is linear between 0.3 
and 1.8. The angular size in all images is 2 arcmin corresponding 
approximately at 10 kpc at the source redshift (z=0.004). 



4.1 Detectability of "outliers" 

As already mention, the heavier the BH with respect to its 
host, the stronger the thermodynamical perturbations are. 
We recall that "heavy" in this context means that the BH 
is above the Mbh — Mbuige relation (i.e. an "outlier"). 

A heavier BH perturbs the background more strongly. 
Figure 7 shows photon images of the X-ray emission de- 
tectable by Chandra, for run D. In this case the presence of 
the BH results not only as a clear distortion of the iso-photal 
contours, but also as a peaked emission separate from the 
centre. This spherical X-ray enhancement survives for many 
passages through the centre. In Figure 8 we clearly distin- 
guish the strong feature of the Mach cone for run E. The 
mirror images, in this case, present sharp well defined edges 
limiting the Mach cone which reaches a maximum contrast 
of 300% over the background. The first panel shows the X- 
ray images when the cone axis is aligned with the horizontal 
axis. It is important to asses if the features are visible along 
lines of sight non exactly orthogonal to the BH motion. We 
produce maps where the angle between the cone axis and 
the plane orthogonal to the line of sight increases from 
(the reference value) to 60 degrees in steps of 5 degrees. The 
signature is recognisable as a cone when the angle is below 
45 degrees, in good agreement with the result we found for 
the density 2D maps alone. For angle above 45 degrees we 
still detect a strong inhomogeneity but the cone structure is 
lost. 

Chandra can resolve the wake of "outliers" to a greater 
distance and in smaller galaxies. The Mach cone is not only 
a key feature for revealing a moving BH with high recoil 
velocities, but it also highlights the presence of an outlier. In 
the following, we rescale our results to galaxies with smaller 
gas fraction and to smaller galaxies, in both cases hosting 
massive BHs above the Mbh — Mbuige relation, "outliers" . 



We first consider how the observability of the conical 
wake depends on the gas fraction. As discussed above, the 
dynamical evolution of the BH is mainly driven by the in- 
teraction with stars, so that the shape of the X-ray features 
is a function only of the properties of the recoiling BH. In 
galaxies with a smaller gas fraction it is only the total X-ray 
flux that changes. 

In order to determine if the wake is visible relative to 
the background galaxy, the corresponding signal to noise^ 
ratio (S/N) needs to be greater than a given detectability 
threshold. The counts expected from a Chandra observation 
of 10 ks give values of the S/N ratio greater than 10 for gas 
fractions above /gas = 10~^, reassuring that, at the distance 
of Virgo, ^ 10^ Mq BHs ejected with supersonic velocity are 
all detectable. 

We now consider the sensitivity of our results to the 
size of the galaxy. Our simulations are scale-free, as long as 
the units we use are consistent with the same gravitational 
constant adopted in GADGET to evolve the system. If we 
scale the masses ( M' = aM), the radii need to be rescaled 
(R' = a'^'^^R, see Shen et al. 2003) in order to preserve 
the relationship between the scale radius and the mass of 
the galaxy: accordingly the time follows t' — a^'^^t. For the 
new system the luminosity of the bremsstrahlung emission 
is related to the old one by L' = a°"^^L. 

In order to determine the minimum value of a for which 
we can still observe the wake, we have to take into account 
not only the S/N ratio, but also the angular resolution of 
Chandra. For our reference galaxy the maximum radius of 
the wake in run b at the first pericentre is ~ 2 kpc. We 
consider the wake in a galaxy at a distance d as resolved 
by Chandra if the angular extent of the wake is covered by 
4 resolution elements (0.5 arcsec, i.e. the angular resolution 
of the telescope). For galaxies at distances comparable the 
Virgo cluster the critical a for a wake to be resolved is aum ^ 
8 X 10~^, corresponding to a BH mass greater than 5 x 
10^ M0. Fixing we have carried on the same analysis 
outlined above for /gas =0.1, 10~^ and 10~^. For /gas =0.1 
and 10~^ the ratio of S/N is again greater than 10. For 
/gas = 10 ~^ it drops to 2 only. 

What happens if we consider more distant galaxies? The 
resolution of Chandra would allow to resolve the wake in 
our reference galaxy up to a distance of ^ 300 Mpc. The 
minimum value of a is related to the distance dhy a ^ 8 x 
10~^((i/20Mpc)^/°-^^. The detection limits imposed by the 
S/N ratio are less tight than those given by the requirement 
on the spatial resolution unless /gas < 10~^. 



5 SUMMARY AND DISCUSSION 

In this paper we studied the perturbation that a recoiling 
massive BH can imprint on the hot gaseous component of 
its host elliptical. The BH is able to affect not only the 
stellar component (as previously noted by Boylan-Kolchin 
et al. (2004) and Gualandris & Merritt (2008)), but also 
the gaseous one. Indeed the ejection of the central object 
perturbs the gas both in its equilibrium properties and along 
its trail. 

^ Note that in this case the noise is given by the emission of the 
host galaxy itself. 
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The BH triggers the formation of a density core on the 
scale of a few rinf, resulting from its sudden removal at the 
time of its ejection and from the kinetic energy deposition 
that continue eroding the central density over the whole 
orbital decay. Along its path, the BH imprints in the gas 
peculiar features that help to reveal its presence while trav- 
elling across its host. The strength of these features depends 
not only by the ejection velocity, but also on the BH mass 
and thus on the relationship between a hole and its host. 
When the motion is subsonic, the BH is surrounded by a 
nearly spherical distribution of gas particles, resulting in an 
enhancement of the underlying bremsstrahlung X-ray emis- 
sion. For the lighter BH (on the Mbh — Mbuige relation), 
the X-ray feature is visible out to the distance of Virgo. 
When the motion is supersonic, the BH shocks the gas ex- 
citing a conical density and temperature perturbation (the 
Mach cone), but the enhanchment does not lead to a de- 
tectable X-ray signal. By contrast, an exceptionally massive 
BH (above the Mbh — Mbuige relation) carries features that 
are not only more intense and detectable in a larger cos- 
mic volume, but that also have a unique flavour: the Mach 
cone is clearly recognisable in X-ray maps with an increase 
in surface brightness by a factor 2-3 over the background. 
Therefore, detection of the Mach cone becomes of twofold 
importance: i) as a probe of high- velocity recoils and ii) as 
an assessment of the scatter of the Mbh — Mbuige relation 
at large BH masses. The different shapes in the disturbance 
can thus in principle help in constraining the extent of the 
gravitational wave recoil and the scatter in scale relationship 
between BH masses and their hosts. 

The lifetime of the different features correlates with 
their shape. For initially low Mach numbers, the spherical 
over-density surrounding the BH, appears every apocentres 
(i.e. where the BH resides most) where the density contrast 
is the highest. This signature lasts up to a hundred Myr. For 
initially high Mach numbers, the Mach cone, recognisable 
during supersonic passages through the centre, is shorter 
lived. Despite its longer duration, the subsonic modes car- 
ries less information as it can be created either by an ini- 
tially fastly moving BH once its orbit has decayed, or by an 
initially slow BH. 

It is important to notice that the X-ray features recog- 
nised in this work are not due to the direct emission from 
the BH but only to the perturbed large-scale hot gas of the 
underlying galaxy. In addition to this, the BH itself can pro- 
duce it own optical and X-ray emission if it is able to retain a 
punctured accretion disc at the moment of its ejection from 
the centre, as suggested by Loeb (2007) (and studied in a 
statistical context by Volonteri & Madau 2008) . The kicked 
BH in this case can turn on as a displaced QSO/AGN, and 
a key manifestation of a large recoil should be imprinted in 
the line systems shifted in velocity from the host galaxies 
(Bonning et al. 2007). This off- nuclear QSO emission can 
coexist with the large scale emission features studied in this 
paper. The contemporary observation of a set of off-set lines 
and a Mach cone can potentially provide information on the 
tangential and radial velocities resulting from the kick. Off- 
centered flares occurring from tidal disruption of bound stars 
or from marginally bound gas that infall on the disc, can also 
be present during the lifetime of the recoiling BH (O'Leary 
& Loeb 2008; Merritt et al. 2008; Shields & Bonning 2008). 

The electromagnetic signature of a recoiling BH, de- 



tailed in this paper, is different from the one expected in the 
immediate vicinity of a coalescence event. Electromagnetic 
afterglows of LISA coalescence events (Milosavljevic & Phin- 
ney 2005; Dotti et al. 2007; Kocsis & Loeb 2008) have pecu- 
liar off-on time depentent features varying on much shorter 
times and spatial scales. 

The detect ability of X-ray features favours heavier BHs 
as the strength of the perturbation increases with Mbh- Ac- 
cordingly, the observability is strongly biased toward more 
massive ellipticals. Interestingly, these massive BHs are out 
of the LISA frequency sensitivity, so that this study would 
complement our knowledge of gravitational wave in-spiral 
events in a window not accessible to LISA. Massive ellipti- 
cals could be the preferred sites where coalescing BHs are 
expected to get the largest kick ever, due to the lack of a 
mechanism able to align the two spins prior to the merger 
(Bogdanovic et al. 2007). The discovery of a prominent X- 
ray feature in the form of a Mach cone will in principle 
support the idea that high velocity recoils are allowed in en- 
vironments where a cool dissipative gaseous component is 
absent. 
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